List of Contents 



Commands 

8/21/72 
4/18/72 
7/11/74 
1/27/72 
2/5/73 

Subroutines 

10/30/72 
9/20/71 



Multics Programmers' Manual 
Installation-Maintained Library 
7/11/74 

bmd 

convert_J3 6 0__f ortran 
list__help 
online_consultant 
troub le_repor t 

f orssp 
plssp 



MULT ICS PROGRAMMERS 1 MANUAL 




Command 

Installation-Maintained Library 

8/21/72 



Name : bmd 



Twelve of the BMD (BioMeDical) programs are available on 
Mul t ics . They are statistical routines developed by the UCLA 
Health Sciences Computing Facility and written in Fortran IV. 
Complete instructions for usage of the programs are available in 
"Biomedical Computer Programs/ 1 edited by W. J. Dixon, University 
of California Publications in Automatic Computation, No. 2. 
Berkeley: University of California Press, 1971. This memo 
contains a list of the routines available and a brief description 
of each one. 

The programs are in ">1 ibraries>bmd / " and this must be in 
the user's search path (i.e., enter the command 
"ssd >libraries>bmd ,r ). 



Usages 

bmdxxx 

where xxx is the particular routine desired. 



Input for the routine, both control cards and data as 
described in the BMD manual, must be in a file called "file07". 
The control cards must be in the fixed format the manual 
specifies. Keywords in control cards, such as "PROBLM, 11 must be 
in all lower case letters, rather than upper case as the manual 
describes. The user may wish to keep the control cards in a 
separate file from the data. In this case, merge the data into 
the proper place in the control file using the merge subcommand 
of edm before execution. Errors in input of this file by the BMD 
programs may result if the file is not closed, so the user must 
be sure the file is properly closed before execution of any BMD 
program. 
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Several of the routines create scratch files as they run. 
The user should be careful not to have files with these names. 
After the program is finished these scratch files should be 
deleted. The routines with scratch files and the names of the 
f i les are: 

bmd02d fileOl 
bmd03m file03 
bmd02r f i 1 eOl, f i 1 e02 
bmd06s f i 1 e01 / f i 1 e02 / f i le03 

The following table lists the routines available: 



bmdOld Simple Data Description — computes simple averages and 
measures of dispersion of variables. 

bmd02d Correlation with Transgenerat ion--computes simple 
correlation coeff icients / averages, and measures of 
dispersion on entering or transgenerated variables from 
selected cases. 



bmd03d Correlation with I tern Delet i on—computes correlation 
matrix, ommitting specified values of variables. The 
punch card option is not available. 

bmd07d Description of Strata with Hi stograms—groups data and 
prints histograms for each variable; means, standard 
deviations, and correlation coefficients are computed 
for each group. 

bmd08d Cross-Tabulation with Variable Stacking — computes two 
way frequency tables of integer data. 

bmd03m General Factor Analysis — performs a principal component 
solution and an orthogonal rotation of the factor 
matr i x. 

bmd05m Discriminant Analysis for Several Groups — computes a 
set of linear functions to classify cases into groups. 



bmd0 2r 



Stepwise Regression — computes a 
linear regression equations in a 



sequence of multiple 
stepwise manner. 
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bmd06s This program performs Guttman Scale analysis. bmd07s 
and bmd08s have been made subroutines of bmd06s. The 
punch card option is not available. The first line in 
the data file should be as follows: 

Col, 1 7 if the program should continue with 
bmd07s / otherwise 0 to stop after bmd06s. 

Col. 2 8 if the program should continue with 
bmd08s, otherwise 0 to stop after bmd07s. 

The next line of data should be the PROBLM card as 
described in the BMD manual. 

bmd02t Autocovar iance and Power Spectral Anal ysi s--computes 

autocovar i ance, power spectrum / cross-covar iance/ 

cross-spectrum, transfer function, and coherence 
function of time series. 
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Command 

Installation-Maintained Library 

4/26/72 



Name: convert 360 fortran, c360f 



This command converts an IBM FORTRAN IV program as read in 
by the Multics card reader to a form acceptable to the Multics 
FORTRAN compiler. 

Usage : 

convert_360__fortran pathl path2 

Pathl is the path or entry name of the segment to be converted. 
This segment name must end in the suffix ".FORTRAN", but the suffix 
need not be typed in the command. path2 is optional and, if 
supplied, is the name of the resulting converted segment. This 
name must have the suffix ".fortran", but, again, the suffix need 
not be typed. If path2 is omitted, the output segment will have 
the same name as the input, except that the suffix will be changed 
from ".FORTRAN" to ".fortran". 

Operation : 

convert_360_fortran makes the following changes to the input 
segment to produce the output segment: 



1) All upper case letters appearing outside of literals of com- 
mentary are changed to lower case. 

2) All single quotes (') appearing outside of literals are changed 
to double quotes ("). Pairs of single quotes appearing in 
literals delimited by single quotes are changed to one single 
quote; double quotes appearing in such literals are replicated. 
Single quotes appearing in literals delimited by nH are not 
changed . 

3) Trailing blanks and the contents of what were card columns 
73-80 are eliminated. Trailing blanks which are a part of a 
literal are preserved. 



4) The card continuation convention is changed to the Multics 
continuation convention. 



convert 360 fortran 
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5) All declaration statements are moved to a position in the 
segment immediately following the program header, if any. As 
a result, only comments, SUBROUTINE, FUNCTION and BLOCK DATA 
statements will precede them. The statements moved are: INTEGER, 
REAL, DOUBLE PRECISION, COMPLEX, LOGICAL, DIMENSION, and COMMON. 

The effect of this is to ensure that DATA and EQUIVALENCE state- 
ments appear in the position required by MULTICS FORTRAN, that 
is, following the above statements. For purposes or rearrange- 
ment, commentary is assumed to be associated with the next 
following non-comment line. 

6) Length specifications in declarations are removed. REAL* 8 is 
changed to "double precision." 



Restrictions : 



convert __360_fortran makes no claim to translate all the incom- 
patibilities between IBM FORTRAN IV and Multics FORTRAN; it does 
only those translations specifically listed above. Suggestions 
for additional translations will be considered. 



The continuatiion convention translation is effected in most 
cases simply be deleting the continuation character on the input 
card image. When a literal crossing card boundaries is encountered, 
the two card images are combined in a single line so that the 
literal is acceptable to Multics FORTRAN. This is the only case 
where this is done, therefore programs which have variables or 
keywords or operators broken across card boundaries will not be 
translated effectively. 

It should be noted that a program which contains statements 
unnecessarily continued on additional cards may be translated 
incorrectly. Specifically, a statement must be classifiable on 
its first card (e.g., the equals sign of an assignment statement 
must occur on the first card) . 

Also note that the conversion from REAL* 8 is not rigorous. 

REAL* 8 XX, YY becomes "double precision xx, yy" ; REAL FUNCTION 

ZZ*8 becomes "double precision function zz" ; but REAL XX, YY*8, 
ZZ*4 becomes "real xx, yy, zz." 
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Command 

Installation-Maintained Library 

07/11/74 



Names : list help, lh 



The list__help command may be used to find help (info) files 
on specified topics. 

Usage : list_help argl . . . argn 

All arguments which are not control arguments are taken as 
"topics" to be searched for. All files whose names contain one 
of the topics will be listed. 

Control arguments : 

-all list the names of all help files 

rather than looking for topics 

-brief list primary names only 

-pathname xxx search specified directory instead of the 

directories used by the help command 

Examples : 

To print the names of files on fortran: list__help fortran 
To look for lisp and macsyma : lh lisp macsyma 

To list all info files in >udd>A>B: lh -pn >udd>A>B -al 

Notes : 

1) The default directories to be searched are: 

>documentation>info segments 
>documentation>iml_J[nf o__segments 

2) The names of help (info) files always end with ".info", but 
list__help does not print this suffix, nor is it specified wh 
printing the contents of a file using the help command. 
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Command 

Installation-Maintained Library 

1/27/72 



Name: online consultant, olc 



This experimental command allows a user to type a message at 
his terminal which will be transmitted to a logged-in IPC staff 
member immediately. The staff member may then reply to the user 
via the "ipc_message_facility . " 

An IPC staff member is scheduled to be logged in every weekday 
between the hours of 2 p.m. and 5 p.m. This schedule will be 
modified if user demand and personnel availability dictate. Out- 
side of the scheduled hours, staff members may sign on as consultant 
(they are encouraged to <io so whenever convenient) , but the Center 
cannot promise an on-line consultant outside scheduled hours. 

Usage : 

There are two ways to use this command: 
online consultant argument (s) 



or 

online_consultant 

If arguments are supplied, a single message consisting of the 
arguments, concatenated and separated by blanks, will be transmitted 
to the current consultant. 



If no argument is supplied, online_consultant will respond 
Input 



and read lines from the console. Each line will be sent to the 
current consultant as a separate message as soon as it is completed. 
The consultant's replies may, of course, be interspersed between 
message input lines. To exit from this message-input mode, the user 
types a line consisting of only a period {.) . 

Notes : 



The on-line consultant will be attempting to complete his normal 
Multics work while signed on. Sometimes, this may mean that the 
messages sent to him are delayed until his program completes execu- 
tion, so please be patient if your question is not answered right 
away . 



online consultant 
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If no consultant is signed on, the command will determine 
whether a consultant is scheduled or not. If not, a short message 
describing other means of bug reporting will be typed. If a 
consultant should be signed on, the command will attempt to locate 
a backup consultant or will type a message giving a telephone 
extension to call. 



The online_consultant command calls the "ipc_jnessage_facility" 
command " accept__messages" as it is exiting, to allow the consultant 
to reply. See the AML writeup for " ipc_message_f acility" for details. 
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Command 

Installation-Maintained Library 

2/5/73 



Names : trouble_report , tr 

The command trouble_report is used to report usage problems 
to the programming consultants at the Information Processing 
Center. Invoking the command provides a series of prompting 
messages designed to assure that all the necessary information 
is reported. The brief option, with only one prompting message, 
is for experienced users of trouble_report. 

Usage ; trouble_report or tr 

This will cause pertinent questions concerning the user and 
the problem to be typed out on the console. The user should 
respond with the requested information or answer yes or no and 
then hit the carriage return once. If the user fails to give 
the requested information, he will receive additional prompting. 

After trouble_report asks the user to describe his programming 
problem, it accepts lines from the console until only a period (" .") 
is typed at the beginning of a line. The user is requested to 
send any source and assembly listings or debugging attempts to 
Consultant, Room 39-454. The path names of the problem segments 
should be given and the ACL should be set to re * . * . * 

Usage : tr -brief or tr -bf 

This types only "Input." Tr then accepts lines from the con- 
sole until a period {" .") is typed at the beginning of a line. 

The temporary segment containing all the positive responses 
to the questions is sent to Consultant IPC by using the mail command. 
The segment is deleted. 

If deemed expedient, the reply to the trouble report will be 
sent to the user's mailbox. In order to permit this, the user 
should see that he has a segment in his login directory named 
"mailbox" and check that the ACL is set to rewa *.*.* If a user 
wants to restrict the access to his mailbox, the following projects 
should have at least "rwa" access to the mailbox and at least "e" 
access to the "home" directory: SysMaint, Multics , PDO, SysLib, 
and Consult. 

When a user wants to restrict the access to the path names of 
segments that can be used in debugging the problem, he should pro- 
vide "re" access to all the project groups mentioned in the previous 
paragraph on the ACL of the directory and on each relevant segment. 
As an alternative, he could set the CACL to "re" if many segments 
are involved. 
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10/30/72 



Name : forssp 



Sixteen of the FORTRAN Scientific Subroutine Package (SSP) 
groups of subroutines are available on Multics. They are Data 
Screening, Correlation and Regression, Design Analysis, 
Discriminant Analysis, Factor Analysis, Time Series, 
Nonparametr i c Statistics, Generation of Random Variates-- 
Distribution Functions, Elementary Statistics and Miscellany, 
Matrices: Storage, Matrices: Operations, Matrices: Inversion, 
Systems of Linear Equations and Related Topics, Matrices: 
E i genana 1 ys i s and Related Topics, Polynomials: Operations, 
Polynomials: Roots, and Polynomials: Special Types. 

Instructions for usage of the programs are available in the IBM 
Svs tern/360 Sc i en t i f i c Subrout i ne Package Programmer 1 s Manua 1 
(GH20-0205). This memo provides a listing of all the subroutines 
(in the same order as in the SSP manual) with a brief description 
of the function of each one. 

Double precision versions have been created for most of the 
routines which did not already have them. These routines are 
starred in this memo and are not found in the SSP manual. They 
work exactly like their single precision counterparts except for 
the typing of their arguments. All real arguments of these 
routines must be double precision except as indicated in the 
following descriptions. 

The subroutines are in ">1 i brar i es>forssp", and this 
must be in the user's search path (i.e., enter the command 
"ssd > 1 i brar i es>forssp") . 



Usage: 

call subrou ti ne_name (arg) 
The subroutines available are: 



Data Screening 



tal ly 
absn t 
submx 



*dtal ly 
*dabsnt 
*dsubmx 



bound *dbound subst *dsubst 
tabl *dtabl tab2 *dtab2 
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cor re 
mu 1 tr 
probt 



Correlation and Regression 



*dcorre 
*dmul tr 
♦dprobt 



mi sr 

gdata 

canor 



*dmi sr 
*dgda ta 
*dcanor 



order 
stprg 



*dorder 
*dstprg 



avdat 



*davdat 



Design Analysis 
avcal *davcal 



meanq 



*dmeanq 



dmatx 



Discriminant Analysis 
*ddmatx disc *ddisc 



trace 



*dtrace 



Factor Analysis 
load *dtoad 



varmx 



*dvarmx 



auto 
exsmo 



*dauto 
*dexsmo 



Time Series 
cross *dcross 



smo 



*dsmo 



Non-Parametr i c Statistics 



kolmo 

krank 

rank 

tie 

wtest 



kolm2 
♦dkrank 
*drank 
*dtie 
♦dwtes t 



smi rn 
mpa i r 
s i gnt 
twoav 



*dsmi rn 
*dmpa i r 
*ds ignt 
*d twoav 



ch i sq 
qtest 
srank 
utes t 



*dchi sq 
*dqtest 
*dsrank 
*dutes t 



General Random Variates — Distribution Functions 



gauss 
*dndtri 



ndtr 



*dndtr 



bdtr 



cdtr 



ndtr i 



mo men 
ph i 
srate 



Elementary Statistics and Miscellany 



*dmomen 
*dph i 
♦dsrate 



ttest 
po i n t 



*dttest 
*dpoi nt 



bi ser 
te tra 



*db i ser 
*dtetra 
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Matrices: Storage 



mcpy 


*dmcpy 


rcpy 


*drcpy 


ccpy 


*dccpy 


dcpy 


*ddcpy 


xepy 


*dxcpy 


ms tr 


*dms tr 


1 oc 


convt 


array 


♦darray 










Matr i ces : 


Opera t ions 






gmadd 


*d.gmadd 


gmsub 


*dgmsub 


gmprd 


*dgmprd 


gmtra 


*dgmtra 


gtprd 


*dgtprd 


madd 


*dmadd 


msub 


*dmsub 


mprd 


*dmprd 


mt ra 


*dmt ra 


tprd 


*dtprd 


ma ta 


*dmata 


sadd 


*dsadd 


ssub 


*dssub 


smpy 


*dsmpy 


sd i v 


*dsd i v 


sc 1 a 


*dsc 1 a 


dc la 


*ddcl a 


radd 


*dradd 


cadd 


*dcadd 


srma 


*dsrma 


scma 


*dscma 


r i nt 


*dr i n t 


ci nt 


*dc int 


rsum 


*drsum 


csum 


*dcsum 


rtab 


*drtab 


ctab 


*dctab 


rsr t 


*drsr t 


csr t 


*dcsr t 


rcut 


*drcu t 


ecu t 


*dccut 


r t i e 


*drt i e 


ct i e 


*dct i e 


mprc 


*dmprc 


mf un 


*dmf un 


recp 


*drecp 



Matrices: Inversion, 
Systems of Linear Equations and Related Topics 



mi nv 


*dmi nv 


s i nv 


ds i nv 


s imq 


*ds imq 


gelg 


dgelg 


rs lmc 


*drs 1 mc 


f actr 


*df actr 


mfgr 


dmf gr 


gel s 


dgel s 


gelb 


dgel b 


mtds 


dmtds 


miss 


dml ss 


mchb 


dmchb 


mf ss 


dmf ss 


mf sd 


dmf sd 


llsq 


dl lsq 




Matr ices : 


E igenanal ys i s and 


Related Topi 


cs 


e i gen 


*de igen 


nroot 


*dnroot 


ateig 


*date i g 


hsbg 


*dhsbg 











Polynomials: Operations 



padd 


*dpadd 


psub 


*dpsub 


pmpy 


*dpmpy 


pdi v 


*dpdi v 


pel a 


*dpcla 


paddm 


*dpaddm 


pva 1 


*dpval 


pvsub 


*dpvsub 


pi Id 


*dpi Id 


pder 


*dpder 


pint 


*dpint 


pqsd 


*dpqsd 


pcld 


*dpcld 


pged 


*dpgcd 


pnorm 


*dpnorm 


pecn 


dpecn 


pecs 


dpecs 
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Polynomials: Roots 



po 1 r t 


*dpol r t 


prqd 


dprqd 


prbm 


dprbm 


pqf b 


dpqf b 












Pol ynomi al s : 


Special 


Types 




cnp 


dcnp 


cnps 


dcnps 


tcnp 


dtcnp 


csp 


dcsp 


csps 


dcsps 


tcsp 


dtcsp 


hep 


dhep 


heps 


dheps 


thep 


dthep 


1 ap 


dl ap 


1 aps 


dl aps 


tlap 


dt lap 


1 ep 


dl ep 


1 eps 


dleps 


tlep 


dtlep 



dlgam 



hist 



Mi seel laneous 
plot 
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Data Screening 



tally to calculate total, mean, standard deviation, 

minimum, and maximum for each variable. 

dtally * double precision version of tally. Only inputs a, 
total, aver, and sd are double precision. 

bound to select the number of observations above, below, 

and between two bounds for each variable. 

dbound * double precision version of bound. 

subst to select a subset from an observation matrix. 

dsubst * double precision version of subst. 

absnt to tabulate missing or zero-values in an observation 

ma tr i x . 

dabsnt * double precision version of absnt. 

tabl to compute frequency tabulation for one variable. 

dtabl * double precision version of tabl. 

tab2 to compute frequency tabulations for two variables. 

dtab2 * double precision version of tab2. 

submx to select a submatrix on the basis of results from 

absnt or subst. 



dsubmx * 



double precision version of submx. 



forssp 
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Correlation and Regression 



corre to calculate means, standard deviations, and 

correlation coefficients. 

dcorre * double precision version of corre. 

misr to calculate means, standard deviations, third 

and fourth moments, correlations, regression 
coefficients, and regression errors. Considers that 
data may be missing. 

dmisr * double precision version of misr. 

order to rearrange i ntercor rel at ions . 

dorder * double precision version of order. 

multr to perform multiple linear regression. 

dmultr * double precision version of multr. 

gdata to generate data matrix for polynomial regression. 

dgdata * double precision version of gdata. 

stprg to perform stepwise multiple linear regression. 

dstprg * double precision version of stprg. 

probt to perform probit analysis. 

dprobt * double precision version of probt. 

canor to perform canonical correlation. 

dcanor * double precision version of canor. 
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Des i gn Anal ys i s 



avdat to prepare data storage for analysis of variance, 

davdat * double precision version of avdat. 

avcal to perform the calculus for the general k-factorial 
experiment: operators 2 and A . 

davcal * double precision version of avcal. 

meanq to perform mean square operation. 

dmeanq * double precision version of meanq. 



Discriminant Analysis 



dmatx to calculate means and dispersion matrix, 

ddmatx * double precision version of dmatx. 

disc to calculate discriminant functions, 

ddisc * double precision version of disc. 



forssp 
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Factor Analysis 



trace to compute cumulative percentage of eigenvalues, 

dtrace * double precision version of trace. 

load to compute factor matrix from eigenvalues and 
associated eigenvectors. 

dload * double precision version of load. 

varmx to perform an orthogonal rotation of a factor matrix, 

dvarmx * double precision version of varmx. 



Time Series 



auto 


to compute autocovar iances. 


dauto * 


double precision version of auto. 


cross 


to compute cross covar iances. 


dcross * 


double precision version of cross. 


smo 


to calculate smoothed or filtered series. 


dsmo * 


double precision version of smo. 


ex smo 


to calculate triple exponential smoothed series. 


dexsmo * 


double precision version of exsmo. 
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Non- Parametr i c Statistics 



kolmo to perform Kolmogorov-Smi rnov one sample test. See 

program comments for double precision usage. 

kolm2 to perform Ko lmogorov-Smi rnov two sample test. See 

program comments for double precision usage. 

smirn to compute values of the Kolmogorov-Smi rnov limiting 

distribution. This program was modified to prevent 
exponent underflows. 

d smirn * double precision version of smirn. 

chisq to compute chi -square for a contingency table. 

dchisq * double precision version of chisq. 

krank to compute the Kendall rank correlation coefficient. 

dkrank * double precision version of krank. 

mpair to perform Wilcoxin's signed ranks test. 

dmpa i r * double precision version of mpair. 

qtest to calculate Cochran's Q. 

dqtest * double precision version of qtest. 

rank to rank a vector of values. 

drank * double precision version of rank. 

signt to perform non-parametric sign test. 

dsignt * double precision version of signt. 

srank to calculate Spearman rank correlation coefficients, 

dsrank * double precision version of srank. 

tie to calculate correction factor due to ties during 
ranki ng. 



forssp 
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dtie * double precision version of tie. 

twoav to calculate Friedman two-way analysis of variance 
stat i st ics . 

dtwoav * double precision version of twoav. 

utest to perform the Mann-Wh i tney U-test. 

dutest * double precision version of utest. 

wtest * to compute the Kendall coefficient of concordance. 

dwtest * double precision version of wtest. 



General Random Variates — Distribution Functions 



gauss to calculate normally distributed random numbers. 

ndtr to compute normal distribution function. 

dndtr * double precision version of ndtr. 

bdtr to compute beta distribution function. 

cdtr to compute chi-square distribution function. 

ndtri to compute inverse of normal distribution function. 

dndtri * double precision version of ndtri. 
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Elementary Statistics and Miscellany 



momen 




to calculate first four moments. 


dmomen 


* 


double precision version of momen. 


ttest 




to test population means. 


dttest 


* 


double precision version of ttest. 


U 1 OC 1 




1"o rnmnnfp h i ^ p r ! A 1 rnr rp 1 Inn rnpf f i r I pn 


dbi ser 


* 


double precision version of biser. 


ph i 




to calculate phi coefficients. 


dph i * 




double precision version of phi. 


point 




to compute po i n t-bi ser i a 1 correlation coefficients 


dpo i n t 


* 


double precision version of point. 


tetra 




to calculate tetrachoric correlation coefficients. 


dte tra 


* 


double precision version of tetra. 


sra te 




to compute survival rates. 


dsrate 


* 


double precision version of srate. 
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Matrices: Storage 



mcpy to copy an entire matrix. 

drncpy * double precision version of mcpy. 

rcpy to copy a specified row of a matrix into a vector, 

drcpy * double precision version of rcpy. 

ccpy to copy a specified column of a matrix into a vector, 

dccpy * double precision version of ccpy. 

dcpy to copy diagonal elements of a matrix into a vector. 

ddcpy * double precision version of dcpy. 

xcpy to copy a portion of a matrix. 

dxcpy * double precision version of xcpy. 

mstr to change the storage mode of a matrix. 

dmstr * double precision version of mstr. 

loc to compute a vector subscript for an element in a 

matrix of specified storage mode, 

convt to convert numbers from single precision to double 

precision or vice versa. 

array to convert a data array from single to double 

dimension or vice versa. This subroutine is used to 
link the user program/ which has double dimension 
arrays, and the SSP subroutines, which operate on 
arrays of data in a vector fashon. 



darray * 



double precision version of array. 
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Matrices: Operations 



gmadd to add two general matrices to form a resultant 

general matrix. 

dgmadd * double precision version of gmadd. 

gmsub to subtract one general matrix from another to form a 

resul tant matrix. 

dgmsub * double precision version of gmsub. 

gmprd to multiply two general matrices to form a resultant 

general matrix. 

dgmprd * double precision version of gmprd. 

gmtra to transpose a general matrix. 

dgmtra * double precision version of gmtra. 

gtprd to premultiply a general matrix by the transpose of 

another general matrix. 

dgtprd * double precision version of gtprd. 

madd to add two matrices element by element to form a 

resul tant matr i x. 

dmadd * double precision version of dmadd. 

msub to subtract two matrices element by element to form a 

resul tant matr i x. 

dmsub * double precision version of msub. 

mprd to multiply two matrices to form a resultant matrix, 

dmprd * double precision version of mprd. 

mtra to transpose a matrix. 



forssp 

Page lk 

dmtra * 
tprd 

dtprd * 
mata 

dmata * 
sadd 

dsadd * 
ssub 

dssub * 
smpy 

dsmpy * 
sdi v 

dsdiv * 
scl a 

dscla * 
dcla 

ddcla * 
radd 

dradd * 
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double precision version of mtra. 

to transpose a matrix and postmultiply by another to 
form a resultant matrix. 

double precision version of tprd. 

to premultiply a matrix by its transpose to form a 
symmetric matrix. 

double precision version of mata. 

to add a scalar to each element of a matrix to form a 
resul tant matrix. 

double precision version of sadd. 

to subtract a scalar from each element of a matrix to 
form a resultant matrix. 

double precision version of ssub. 

to multiply each element of a matrix by a scalar to 
form a resultant matrix. 

double precision version of smpy. 

to divide each element of a matrix by a scalar to 
form a resultant matrix. 

double precision version of sdiv. 

to set each element of a matrix equal to a given 
scalar. 

double precision version of scla. 

to set each diagonal element of a matrix equal to a 
seal ar . 

double precision version of dcla. 

to add a row of one matrix to a row of another 
ma tr i x. 

double precision version of radd. 
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cadd to add a column of one matrix to a column of another 

matrix. 

dcadd * double precision version of cadd. 

srma to multiply a row of a matrix by a scalar and to add 

to another row of the same matrix. 

dsrma * double precision version of srma. 

scma to multiply a column of a matrix by a scalar and to 

add to another column of the same matrix. 

dscma * double precision version of scma. 

rint to interchange two rows of a matrix. 

drint * double precision version of rint. 

cint to interchange two columns of a matrix. 

dcint * double precision version of cint. 

rsum to sum the elements of each row of a matrix to form a 
column vector. 

drsum * double precision version of rsum. 

csum to sum the elements of each column of a matrix to 

form a row vector. 

dcsum * double precision version of csum. 

rtab to tabulate rows of a matrix to form a summary 

matr i x. 

drtab * double precision version of rtab. 

ctab to tabulate columns of a matrix to form a summary 

matrix. 

dctab. * double precision version of ctab. 

rsrt to sort rows of a matrix. 

drsrt * double precision version of rsrt. 
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csrt 
dcsrt * 
rcut 

drcut * 
ccut 

dccut * 
r t i e 

drtie * 
ct i e 

dctie * 
mprc, dmprc 

mf un 

dmfun * 
recp 



to sort columns of a matrix, 
double precision version of csrt. 

to partition a matrix between specified rows to form 
two resultant matrices. 

double precision version of rcut. 

to partition a matrix between specified columns to 
form two resultant matrices. 

double precision version of ccut. 

to adjoin two matrices with the same column dimension 
to form one resultant matrix. 

double precision verison of rtie. 

to adjoin two matrices with the same row dimension to 
form one resultant matrix. 

double precision version of ctie. 

to permute the rows or columns of a given matrix 
according to a given transposition vector or its 
i nverse . 

to apply a function to each element of a matrix to 
form a resultant matrix. 

double precision version of mfun. 

to calculate the reciprocal of an element. This is a 
FORTRAN function subprogram that may be used as an 
argument by the subroutine MFUN. 



drecp * 



double precision version of recp. 
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Matrices: Inversion, 
Systems of Linear Equations and Related Topics 



mi nv 
dm i n v * 
s i nv, dsl nv 
s imq 
dsimq * 
gelg,dgelg 

rs 1 mc 

drslmc * 

f act r 

dfactr * 
mf gr, dmf gr 
gel s, dgel s 

gel b, dgel b 

mtds, dmtds 
ml ss, dml ss 

mchb, dmchb 



to invert a square matrix. 

double precision version of minv. 

to invert a symmetric positive definite matrix. 

to solve simultaneous linear algebraic equations. 

double precision version of simq. 

to solve a system of general simultaneous linear 
equations by Gauss elimination. 

to solve simultaneous linear equations with iterative 
ref i nement . 

double precision version of rslmc. Inputs a, af, b, 
x, v, and per should be double precision. Input epsi 
should be single precision. 

to compute a triangular factorization of a 
nonsingular matrix. 

double precision version of factr. 

to factor a matrix and determine its rank. 

to solve a system of general simultaneous linear 
equations with symmetric coefficients. 

to solve a system of general simultaneous linear 
equations with band-structured coefficients. 

to divide a matrix by a triangular matrix. 

to solve simultaneous linear equations with symmetric 
positive semidefinite matrix. 

to compute a triangular factorization of a symmetric 
positive definite band matrix. 
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mfss,dmfss to compute a triangular factorization and determine 
the rank of a symmetric positive semidefinite matrix. 

mfsd/dmfsd to compute a triangular factorization of a symmetric 
positive definite matrix. 

llsq,dllsq to solve linear least squares problem. 



Matrices: E i genana 1 ys i s and Related Topics 



eigen to compute eigenvalues and eigenvectors for a real 

symmetric matrix. 

deigen * double precision version of eigen. 

nroot to calculate eigenvalues and eigenvectors for a real 

nonsymmetric matrix of the form B-inverse times A. 

dnroot * double precision version of nroot. 

hsbg to reduce a real matrix into upper almost triangular 

form (hessenberg form). 

dhsbg * double precision version of hsbg. 

ateig to compute eigenvalues for a real almost triangular 

ma t r i x . 



dateig * 



double precision version of ateig. 
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Polynomials: Operations 



padd to add two polynomials, 

dpadd * double precision version of padd. 

psub to subtract two polynomials, 

dpsub * double precision version of psub. 

pmpy to multiply two polynomials, 

dpmpy * double precision version of pmpy. 

pdiv to divide one polynomial by another. 

dpdiv * double precision version of pdiv. Inputs p, x, and y 

should be double precision. Input tol should be 
single precision. 

pcla to replace one polynomial by another. 

dpcla * double precision version of pcla. 

paddm to multiply one polynomial by a constant and add it 

to another polynomial. 

dpaddm * double precision version of paddm. 

pval to evaluate a polynomial. 

dpval * double precision verision of pval. 

pvsub to substitute the variable of one polynomial by 

another polynomial. 

dpvsub * double precision version of pvsub. 

pild to evaluate a polynomial and its first derivative. 

dpild * double precision version of pild. 



forssp 
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pder 
dpder * 
pint 
dpint * 
pqsd 

dpqsd * 
pcld 

dpc 1 d * 
pgcd 

dpgcd * 

pnorm 
dpnorm * 

pecn / dpecn 
pecs, dpecs 



to take the derivative of a polynomial, 
double precision version of pder. 
to integrate a polynomial, 
double precision version of pint. 

to perform quadratic synthetic division of a 
pol ynomi al . 

double precision version of pqsd. 

to perform complete linear synthetic division of a 
pol ynomi al . 

double precision version of pcld. 

common 



to find the greatest 
pol ynomi al s . 



double precision version of pgcd. 
work should be double precision, 
single precision. v 



divisor of 



two 



Inputs x, y, and 
Input eps should be 



to normalize the coefficient vector of a polynomial. 

double precision version of pnorm. Only input x 
should be double precision. Input eps should be 
single precision. 

to economize a polynomial for a symmetric range, 
to economize a polynomial for an unsymmetric range. 
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Polynomials: Roots 



polrt to find real and complex roots of a real polynomial, 

dpolrt * double precision version of polrt, 

prqd / dprqd to find roots of a real polynomial by QD algorithm 
with displacement. 

prbm / dprbm to find roots of a real polynomial by Bairstow's 
a 1 gor i thm. 

pqfb / dpqfb to determine a quadratic factor of a real polynomial. 
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Polynomials: Special Types 



cnp,dcnp to compute the values of the Chebyshev polynomials 
for a given argument x and orders zero up to n. 

cnps,dcnps to compute the value of a series expansion in 
Chebyshev polynomials. 

tcnp,dtcnp to transform a given series expansion in Chebyshev 
polynomials to a polynomial. 

csp/dcsp to compute the values of the shifted Chebyshev 
polynomials for a given argument x and orders zero up 
to n. 

csps,dcsps to compute the value of a series expansion in shifted 
Chebyshev polynomials. 

tcsp / dtcsp to transform a given series expansion in shifted 
Chebyshev polynomials to a polynomial. 

hep/dhep to compute the values of the Hermite polynomials for 
a given argument x and orders zero up to n. 

hepS/dheps to compute the value of a series expansion in Hermite 
pol ynomi als. 

thep/dthep to transform a given series expansion in Hermite 
polynomials to a polynomial. 

lap/dlap to compute the values of the normalized Laguerre 
polynomials for a given argument x and orders zero up 
to n. 

laps,dlaps to compute the value of a series expansion in 
normalized Laguerre polynomials. 

tlap/dtlap to transform a given series expansion in normalized 

Laguerre polynomials to a polynomial. 

lep,dlep to compute the values of the Legendre polynomials for 
a given argument x and orders zero up to n. 
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leps,dleps to compute the value of a series expansion in 
Legendre polynomials. 

tlep/dtlep to transform a given series expansion in Legendre 
polynomials to a polynomial. 



Mi seel laneous 



dlgam to compute the double precision natural logarithm of 

the gamma function of a given double precision 
argument. This routine is from the Special Functions 
section and is used by routines bdtr and edtr from 
the Random Variates section. 

hist to print a histogram of frequencies versus intervals. 

This routine is from the Special Sample Subroutines 
sect ion . 

plot to plot several cross-variables versus a base 

variable. This routine is from the Special Sample 
Subroutines section. 
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Subroutine 

Installation-Maintained Library 

9/20/71 



Name ; plssp 

All but 15 members* of the PL/1 Scientific Subroutine Package 
(SSP) are available on Multics. Instructions for using the sub- 
routines are available in the IBM System/360 Scientific Subroutine 
Package (PL/1) Program Description and Manual (GH20-0586-0) . 

This memo contains directions for using the subroutines on 
Multics and a list of all the subroutines available with a brief 
description of the function of each one. 

The subroutines are in ">libraries>plssp" , and this must be in 
the user's search path (i.e., type "ssd >libraries>plssp") . 

Usage : 

call subroutine_name (arg) 

The subroutines available and descriptions of the functions of 
these subroutines are on the following pages. 



*ahie ,ahim,bdtr ,cdtr , f f t , f f tm, fmfpjmon^mva^mveb/mvs^mvubjpec , 
prtc and ptc will not be available until the installation of the 
new Multics PL/1 compiler. 



plssp 
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1. 



Mathematics 



acfe 

acfm 
alie 

alim 

apcl 

ape 2 

apll 

asn 

cell 

cel2 
dere 



interpolates the function value yval for a given argument 
xval using xst f the starting value of the arguments, dx, 
the increment of the argument values and vector y of func- 
tion values . 

interpolates the function value yval for a given argument 
xval using a given table (x,y) of arguments and function 
values . 

interpolates the function value yval for a given argument 
value xval using xst, the starting value of the arguments, 
dx, the increment of the argument values and the vector 
y of function values. 

interpolates the function value yval for a given argument 
value xval using a given table (x,y) of argument and function 
values . 

sets up the normal equations for a weighted polynomial least 
squares fit to a given discrete function, using Chebyshev 
polynomials as fundamental functions . 

sets up the normal equations for a polynomial least squares 
fit to a given discrete function, using Chebyshev poly- 
nomials as fundamental functions. 

sets up normal equations for a polynomial least squares fit 
to a given discrete function. 

computes the solution of normal equations set up by pro- 
cedures apcl, apc2 and apll. 

computes the complete elliptic integral of the first kind: 



computes the generalized elliptic integral of the second 
kind: 



performs one integration step for a system of first order 
ordinary differential equations y'=f(x,y) with given initial 
values y. The stepsize h is adjusted for accuracy require- 
ments and speed considerations. 
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det3 computes a vector z = (z 2 , z £ , . .., z^^) of derivative 

values, given a vector y = (y l , y 2 , — , y dim ) of function 

values whose components y^ correspond to dim equidistantly 

spaced argument values x. with x. - x. , = h for i = 2, 
dim. 1 1 1 ~ 1 

det5 computes a vector z = (z x , z 2 , z^ m ) of derivative 

values y = (y x , y 2 , y^j^) of function values whose 

components correspond to dim equidistantly spaced argument 
values x^, with x^ - x^_ 1 = h. 

dfeo given argument x and function fct(x) , defined in the one- 
sided interval [x,x+h] , dfeo computes an approximation 
z to the derivative. 

dfec given the argument x and the function fct(x) , defined in 
the closed interval [x-|h|, x+|h|], dfec computes an 
approximation z to the derivative of the function fct(x) . 

dgt3 computes a vector z - (z jf z 2 , z ^ m ) / °f derivative 

values, when vectors x = (x 1# x 2 , x^j^) of argument 

values and y = (y lf . y<ji m ^ °f corresponding 

function values are given. 

elil computes the incomplete elliptic integral of first kind for 
given values of an argument x and complementary modulus ck. 

eli2 computes the incomplete elliptic integral of second kind 

for given values of an argument x, complementary modulus ck 
and constants a and b. 

exsm develops the triple exponential smoothed series s of the 
given series x. 

jelf calculates the three Jacobian elliptic functions sn, cn and 
dn. 

lgam computes the double-precision natural logarithm of the 
gamma function of a given double-precision argument. 

post transforms a given series expansion in orthogonal polynomials 
to a polynomial. The independent variable of the given 
expansion is assumed to be x 0 + x^x; that is, a linear 
trans foramtion of the range is built in. The coefficient 
vector c - (c^ c^, c ) is given. Procedure post 

calculates pol = (pol^ pol 2 , P°l n ) satisfying 

n n 



^ C i f i-l (x o +X i* x)= P ( P° 1 i ) (xi ~ 1) 
i=l i=L 
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For the specified set of orthogonal polynomials (f^) the 
user has the choice of : 

Chebyshev polynomials (t Q , t x , . .., t ^) with opt = "t" 

Legendre polynomials (p 0 , p x , p with opt = "p" 

Laguerre polynomials (1 0# l jr l n-1 ) with opt = "1" 

Hermite polynomials (h Q , h Jf h n _]_) with opt = "h" . 

computes the value of the sum 
n 

c k f k-l (x) 
k=l 

for a given vector c = (c l# c 2 , c n ) and a specified 

set of orthogonal polynomials (f^) . 

The user has a choice of: 

Chebyshev polynomials (t 0 , t jr . .., t ^) with opt = "t" 

Legendre polynomials (p Q , p 1# . .., P n .^) with opt = "p" 

Laguerre polynomials (1 Q , 1 x , . .., ^ n «i^ with opt = "1" 

Hermite polymonials (h Q , h jf . .., n n »i^ with opt = "h" . 

computes the values of the first n orthogonal polynomials . 
The user has the choice of ; 

Chebyshev polynomials (t Q , t Jf t ^) with opt = "t" 

Legendre polynomials (p flf p lf P n _^) with opt = "p" 

Laguerre polynomials (1 0 , l jf . .., ^ n -i^ with opt = "1" 

Hermite polynomials (h fl , h lt . .., n n _i) with opt = "h" . 
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se13 computes a vector z = (z t , z 2 , z d ^ m ) of smoothed func- 

tion values. It requires a vector y = (y x , y 2 , y^^ m ) 

where the y components correspond to equidistantly spaced 
argument values x^, assuming x^ - - h. 

se15 computes a vector z = {z x , z 2 , . .., z ^i m ^ °f smoothed func- 
tion values, given a vector y = (y jf y 2 , . .., Y^im^ °f func- 
tion values whose components y. correspond to dim equi- 

x 

distantly spaced argument values x^ with x^ - = h for 

i = 2 , . . . , dim. 

se35 computes a vector z = (z jf z , . .., z^^ m ) of smoothed 
function values, given a vector y = (y t , y 2 , Y^^) 
of function values whose components y^ correspond to 
dim equidistantly spaced argument values x^ with 
x^ - x ^_i = n f° r i = 2, dim. 



sg13 computes a vector z = (z x , z 2 , z d ^ m ) ° f smoothed 

function values. It requires a vector y = (yj , y 2 , 
and a vector x = (x Jf x 2 , x dim ) of argument values, 

where y^ corresponds to x^. 



dim' 
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2. Matrix Operations 

cane computes the canonical correlations between two sets of 
variables. It is normally preceded by a call to pro- 
cedure corr. 

mags computes c = a + b if opt = "1" 

c = a - b if opt = "2" 

c = b - a if opt = "3" 
for given matrices a and b which are general and symmetric 
respectively . 

mate reduces a given real matrix to upper almost triangular 

(Hessenberg) form by means of a sequence of similarities. 

matu reduces a given real matrix to almost triangular (Hessenberg) 
form by means of a sequence of orthogonal transformations. 

mdlg for a system of equations a . x = r, where a = 1 • u is a 
general nonsingular matrix, mdlg performs the following 
calculations , depending on the character of an input parame- 
ter opt: 

opt = "1" r is replaced by l" 1 • r 
opt = "2" r is replaced by u" 1 • r 
otherwise r is replaced by (1 . u)- 1 . r. 

mdls for a system of equations a • x = r with symmetric positive 

definite matrix a = t • t t , mdls performs the following 
calculations depending on the character of the input parame- 
ter opt: 

opt = "1" r is replaced by t" 1 • r 
opt = "2" r is replaced by (t" 1 ) t • r 
otherwise r is replaced by (t • t"*")' 1 • r. 

mdrs for a system of equations x * a - r with symmetric positive 

4- 

definite matrix a = t • t , mdrs performs the following 
calculations , depending on the character of an input parame- 
ter opt: 

opt = "1" r is replaced by r • (t~ 1 ) t 
opt = "2" r is replaced by r • f* 1 

otherwise r is replaced by r • (t • t 1 ")" 1 . 

mdsb depending on the character of the input parameter opt r mdsb 
performs the following operations on a system of equations 
a • x = r with symmetric positive definite band matrix 

a = t fc • t: 

opt = "1" r is replaced by (t"" 1 ) t • r 
opt = "2" r is replaced by t" 1 • r 



otherwise r is replaced by (t • t) -1 
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meat computes the eigenvalues of a real upper almost triangular 
matrix (Hessenberg form, see subroutines mate and matu) 
using double OR iteration. 

mebs computes a lower and an upper bound for the eigenvalues 
of a real symmetric matrix. 

mest computes the eigenvalues of a real symmetric tridiagonal 
matrix (see subroutine mstu) . 

mfg factorizes a general nonsingular matrix a into a product 
of a lower triangular matrix 1 and an upper triangular 
matrix u overwritten on a, omitting the unit diagonal of 
u. 

mfgr performs the following for a given rectangular matrix: 

determines rank and linearly independent rows and 
columns 

factorizes a submatrix of maximal rank 

expresses nonbasic rows in terms of basic rows 

expresses basic variables in terms of free variables. 

mfs computes a triangular factorization of a symmetric positive 
definite matrix using the square root method of Cholesky. 

mfsb computes a triangular factorization of a symmetric positive 
definite band matrix using the square root method of 
Cholesky. 

mgbl performs the following operations on an equation system 
a • x = r with general band matrix a = 1 • u, depending 
on the character of an input parameter opt: 



opt = "1" 


u 


replaces a and 


!-i . 


r 


replaces r 


opt = "u" 


u 


replaces a and 


u- 1 • 


r 


replaces r 


opt = "b" 


u" 


'* • r replaces 


r for 


a 


given u on storage 



locations of a 

otherwise u replaces a and the solution x = a" 1 
replaces r. 
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mgb2 performs the following operations on an equation system 
a • x * r with general band matrix a = 1 • u, depending 
on the character of an input parameter opt: 

opt = "1" a is replaced by upper band factor u, r is 
replaced by l" 1 • r and lower band factor 
1 is stored in a one-dimensional array 1 
omitting the unit diagonal 

opt = "f" a is replaced by the upper band factor u and 
the lower band factor 1 is stored in the 
array 1; the right-hand side r remains 
unchanged 

opt = "a" r is replaced by l" 1 • r for the given upper 
factor u in array a and the lower factor 1 
in vector 1 

opt = "c" r is replaced by the solution x = a~* • r 
for given u and 1 

otherwise a is replaced by the upper factor u; the 
lower factor 1 is calculated and stored in 
1, and r is replaced by the solution x = a -1 • r. 

mgdu computes eigenvalues and eigenvectors of a real matrix 
of the form b-inverse times a, where a is symmetric and 
b is positive definite. 

mig inverts a general nonsingular matrix a, which is given in 
the factored form a = 1 • u, where the upper triangular 
matrix u contains the unit diagonal which is not stored. 

minv inverts a general square matrix. 

mis inverts a symmetric positive definite matrix a, which is 
given in factored form (Cholesky) : 

a = t • transpose (t) . 

mlsq calculates x satisfying a • x = b, that is, the solution 
of a system of linear equations using Householder trans- 
formations. The least squares solution is obtained 
in case of an overdetermined system of equations . 



mmgg 



computes the standard matrix product c = a • b. 
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mmgs calculates g • s if opt = "1" or 

s • g if opt = "2" 
where g is a general and s a symmetric matrix. 
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mmgt calculates a • a if opt = "1" or 

a 1 • a if opt = "2" 
where a is a general matrix. 

mmss computes the standard product p = a 
matrices . 



mpit 



mprm 



b of two symmetric 



Calculates the permutation vector corresponding to the 
inverse of a given permutation if opt = "i" and the trans- 
position vectors of the given permutation and of its 
inverse if opt = "t". 



permutes rows (if opt = "r") or columns (if opt = "c") 

.f 

p * (it mv = "i") . Tne permu- 
tation p is given in the form of its transposition vector 
t. 



of a given matrix a according to the permutation p (i: 
inv = "0") or its inverse p~* (if inv = "1") . The pe; 



mscg expands the compressed one-dimensional storage allocation 
of a symmetric matrix to general two-dimensional form. 

mscs compresses the storage allocation of a symmetric two- 
dimensional matrix to a one-dimensional array. 

msdu computes eigenvalues and eigenvectors of a real symmetric 
matrix . 

mstu reduces a given real symmetric matrix to tridiagonal form 
by means of a sequence of orthogonal transformation. 

mtpi calculates the permutation vector if inv = "0" and the 

inverse permutation vector if inv = "1" from a given trans- 
position vector. 

mvst provides the eigenvector corresponding to a given eigenvalue 
of a symmetric tridiagonal matrix. 
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3. Statistics 

abst tests for missing or zero elements in observation matrix a. 

avar performs an analysis of variance for a complete factorial 
design. 

boun selects from a set (or a subset) of observations the number 
of observations under, between and over two given bounds 
for each variable. 

chsq computes chi-square from a contingency table. 

corr computes means, standard deviations, sums of cross-products 
of deviations and correlation coefficients. 

dmtx computes means of variables in each group and a pooled 

dispersion matrix for all the groups. This subroutine is 
used in the performance of discriminant analysis . 

dscr performs a discriminant analysis by calculating a set of 
linear functions that serves as indices for classifying 
an individual into one of k groups. 

htes calculates the Kruskal-Wallis H-statistic from the ranks 
of observations obtained from three or more independent 
samples . 

klmo tests the difference between empirical and theoretical 
distributions using the Kolmogorov-Smirnov test. 

klm2 tests the difference between two sample distributions using 
the Kolmogorov-Smirnov test. 

krnk measures the correlation between two variables by means 
of the Kendall rank correlation coefficient. 

load computes a factor matrix (loading) from eigenvalues and 
associated eigenvectors. 

mltr performs a multiple linear regression analysis for a depen- 
dent variable and a set of independent variables. 

ndti computes x = p~My) such that y = p(x) , the probability 

that the random variable x, distributed normally (0,1), is 
less than or equal to x; f (x) , the ordinate of the normal 
density at x, is also computed. 
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ndtr computes y = p(x) , the probability that the random avail- 
able x, distributed normally (0,1), is less than or equal 
to x; f (x) , the ordinate of the normal density at x, 
is also computed. 

ordr is used to choose a dependent variable and a set of indepen- 
dent variables from a matrix of correlation coefficients 
and to form a submatrix of correlation coefficients to 
be used in performing a multiple linear regression analysis. 

qtst uses the Cochran Q-test to determine whether three or more 
matched groups of dichotomious data differ sygnif icantly . 

rank ranks a vector of data. 

sbst derives a subset vector indicating which observations in 
a set have satisfied certain conditions on the variables. 

smir computes values of the limiting distribution function for 
the Kolmogorov-Smirnov statistic. 

srnk tests the correlation between two variables by means of 
the Spearman rank correlation coefficient. 

strg performs a setpwise multiple linear regression analysis 

for a dependent variable and a set of independent variables. 

subm copies a submatrix from an observation matrix. The elements 
of this submatrix satisfy conditions specified by an input 
vector. This subroutine is used in preparing data for input 
to a statistical analysis such as multiple regression. 

tab1 tabulates, for one variable in an observation matrix (or a 
matrix subset) , the frequency and percent frequency over 
given class intervals. In addition, t calculates for the 
same variable the total, mean, standard deviation, minimum 
and maximum. 

tab2 performs a two-way classification, for two variables in 

an observation matrix (or a matrix subset) , of the frequency, 
percent frequency and other statistics over given class 
intervals . 

taly calculates total, mean, standard deviation, minimum and maxi- 
mum for each variable in a set (or a subset) of observations. 

tie calculates correction factor due to ties. 



trac 



computes cumulative percentage of eigenvalues greater than 
or equal to a constant specified by the user. 
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ttst calculates certain t-statistics on the means of populations. 

twav tests whether a number of samples are from the same popula- 
tion by the Friedman two-way analysis of variance test. 

utst tests whether two independent groups are from the same 
population by means of the Mann-Whitney U-test. 

vrmx performs an orthogonal rotation of a factor matrix. 

wtst measures the degree of association among a number of variables 
by the Kendall coefficient of concordance. 



